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Chile is expecting a 5.4% growth in energy consumption per year until 2030, requiring new and better 
solutions for the upward trend of its electricity demand. This state leads to select and study one of the 
potential alternatives for electricity generation: concentrated solar power (CSP) plants. Such renewable 
technology found in Chile a very favorable condition. Recent researches indicate Atacama Desert as one 
of the best regions for solar energy worldwide, having an average radiation higher than in places where 
CSP plants are currently implemented, e.g. Spain and USA. The aim of this study is to present an analysis 
of levelized energy cost (LEC) for different power capacities of CSP plants placed in distinct locations in 
northern Chile. The results showed that CSP plants can be implemented in Atacama Desert with LECs 
around 19 ¢US$/kWh when a gas-fired backup and thermal energy storage (TES) systems are fitted. 
This value increases to approximately 28 ¢US$/kWh if there is no backup system. 
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1. Introduction per year until 2030 [1]. The current energy situation in Chile and 


its radiation conditions gives the opportunity to investigate the 


Chile is nowadays passing through a period of uncertainty in 
what regards its energy sector, which is found to be unprepared 
to provide the estimated growth of 5.4% in energy consumption 


Abbreviations: CSP, concentrated solar power; LEC, levelized energy cost; SPT, 
solar parabolic trough; TES, thermal energy storage; NCRE, non-conventional 
renewable energy; SIC, electric central interconnected system; SING, electric 
northern interconnected system; UTFSM, Universidad Técnica Federico Santa 
Maria; NSA, national solarimetric archive; DNI, direct normal irradiance; HTF, heat 
transfer fluid; PB, power block; PTC, parabolic trough collector; DSG, direct steam 
generation; AIC, annual electricity generation; TICVC, total investment costs and 
variable costs 
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potential and viability of implementing solar thermal power 
plants in the north of the country [2-4]. However, there are only 
few researches available in this subject regarding Chile, and the 
majority of them date from two or three decades ago [5-9]. In 
addition to that, it is noted that none of the existing researches 
consider economic aspects to install a concentrated solar plant 
in Chile. 

Regarding Chilean energy regulations, it is important to 
remark that, since 2008, the electricity companies selling directly 
to final customers must incorporate a certain percentage of non- 
conventional renewable energy (NCRE) into the electricity they 
trade [10]. This law consolidates the efforts to remove barriers to 
the incorporation of NCRE into the national power mix, thereby 
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contributes to the objectives of supply security and environmen- 
tal sustainability that govern Chile’s energy policy. Nevertheless, 
according to the latest official report on the subject [10], the 
supply of energy in Chile is primarily provided by conventional 
power plants, as coal-fired thermal plants, combined cycle, diesel 
engines, gas/oil turbines, and large scale hydro electrical power 
plants. As shown in Fig. 1, in 2010, only approximately 3% of the 
energy was produced by non-conventional renewable power 
plants (small hydro, biomass, and wind). In the same year, from 
a total of 14.878 MW of installed capacity (Fig. 2), only 3.1% was 
NCRE plants. 


Biomass: 1% Small Hydro 
aa plants < 
20MW: 1% 


Wind: 1% 


Fig. 1. Electricity generation in Chile (SIC and SING) during 2010 [10]. 
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Fig. 2. Installed capacity in Chile (SIC and SING) by generation technology 
in 2010 [10]. 
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Fig. 3. Forecast of LEC values for different generation technologies [14,13]. 


The statistics show that there is a very favorable scenario for 
absorbing new NCRE projects needed to accomplish the percen- 
tage of NCRE officially required. Therefore, in order to encourage 
private and public investments in solar thermal energy, it is 
necessary for the national academic community to provide novel 
studies and tools in which technical and economic analysis can 
rely on. Concentrated solar power (CSP) technology has been 
widely investigated and enhanced in many countries [11,12], and 
these studies serve as a solid basis for newer investigations of CSP 
technology applied to the Chilean context. Moreover, it has been 
found that, for the next few years, costs of solar thermal energy 
generation present a trend of decrease, while the fossil-fuel 
alternative presents an upward evolution. According to Fig. 3 
[13-15], it is possible for solar thermal energy costs to be at the 
same level of fossil-fuel energy costs around 2020. 

Considering all explained above, the present study aims to 
develop an economic analysis for computing the LEC acquired 
from different CSP plant sizes in distinct locations in northern 
Chile. The research basis lay on contributions of recent Chilean 
researches [2,4] together with international researches [16-21] in 
order to evaluate the implementation of CSP plants in the north of 
the country. The analysis is undertaken through a computational 
simulation, which is performed for two different plant arrange- 
ments: the first consists in a solar parabolic trough (SPT) power 
plant fitted with thermal energy storage (TES) system and a gas- 
fired backup system; and a second configuration of a SPT plant 
equipped only with a TES system. Thereby, the simulation allows 
determining the plant power capacity and configuration that 
provides the lowest LEC among the different locations selected. 


2. Potential for concentrated solar plants in northern Chile 


Recent researches regarding the potential of solar radiation in 
northern Chile [2,4]indicates the Atacama Desert as being one of 
the best worldwide regions for solar energy, based on energy 
density data from several sources [22-25]. According to them, the 
region presents a high number of clear days during the year, due 
to the aridity of the Atacama Desert, defined as a hyperacid region 
[26] with annual average precipitations lower than 50 mm per 
year. This research considered the solar radiation data presented 
by [27] also available at [24], which is derived from ground 
station measurements compiled by Universidad Técnica Federico 
Santa Maria (UTFSM). More discussion about the solar energy 
data in Chile and its uncertainty levels can be found in the 
literature [25]. Table 1 [4,18] and Fig. 4 [27] show that the 
average radiation in the north of Chile is even higher than in 
some places where solar power plants are currently implemented. 

For the economic evaluation of CSP plants, the present study 
accounts not only the Atacama Desert but also the several 
locations in northern Chile within the SING (northern electric 
interconnected system) area. It is considered that different 
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regions in this territorial range have a strong potential for placing 
CSP plants. 

In the following section, the solar radiation data extracted 
from monthly means provided by UTFSM is processed in order to 
create an hourly database. The location and technology selection 
are also presented, as well as the method used for calculating the 
final energy produced from the available solar radiation. Finally, 
an analysis of the Levelized Energy Costs (LEC) is proposed for 
establishing the most suitable plant size and configuration 
required for each region analyzed. 


3. Methodology 
3.1. Estimation of available solar radiation and location selection 


About 16 locations in the Atacama Desert are presented in the 
UTFSM data and chosen for evaluation, all of them situated 
between —18.20° y —26.27° latitudes and —70.46° y —68.01° 
longitudes. 

In the country, four radiation databases are available to the 
public, all of them obtained by different ground station networks: 


e National Solarimetric Archive (NSA): 89 stations. 
e Meteorological service data: 18 stations. 
e CNE-GTZ ground stations: 7 stations. 
e UC-FONDEF ground stations: 12 stations. 
The radiation data used for the simulations were acquired 
from the national solarimetric archive [27], which include data 
covering from 1961 to 1983. In comparison with the additional 
databases, NSA registrations are in between the others, with 
moderately low radiation levels, making it appropriate for eco- 
nomic evaluations in a worst-case scenario. 

In the locations selected, solar radiation peaks during sum- 
mertime at values close to 900 MJ/m? month, with winter minima 


N. Corral et al. / Renewable and Sustainable Energy Reviews 16 (2012) 6678-6685 


nearing 260 MJ/m2 month. Fig. 4 shows global average irradiance 
per month per square meter. 

Chilean solar radiation data are available only in the form of 
monthly means measurements of the daily total radiation (H 
for each location).Therefore, it is necessary to apply a methodol- 
ogy that converts these values into hourly databases, a more 
accurate data to use as an input in the power plant simulation. 
Considering that a CSP plant will only accept direct normal 
irradiance (DNI) in order to work, the main objective is to develop 
a model for calculating the clear days. The procedure used is 
outlined by Duffie and Beckman [23], and is based on the Liu and 
Jordan [28] generalized distributions of cloudy and clear days, 
later modified by Bendt et al. [29,30], then by Stuarts and 
Hollands [31] and finally by Knight [32]. Based on this theoretical 
background, the hourly radiation beam is obtained using a 
correlation method described in [4,23,32,33]: first, an artificial 
daily total radiation H is obtained for each day of the month; then, 
an artificial month of daily and hourly values of solar radiation is 
built; finally, the clear day model is applied to obtain hourly 
values. 


3.2. Estimation of daily radiation 


For converting available data into hourly database, the pre- 
viously mentioned methodology is based in the clearness index 
(Kr), which is the ratio between daily radiation in a horizontal 
surface and the extraterrestrial radiation. The selected procedure 
is also used in [4] and is described essentially in [23]. 

First, it is necessary to compute monthly average clearness 
index (Kr) for each month and location, which is defined as 
= H 
Kr = Ao a) 
where H is the monthly average radiation, obtained from the 
registered measurements in [27], and Hp is the monthly average 
extraterrestrial radiation, computed for each day and location by 
the following formula: 


Table 1 
Solar radiation of CSP plants and projects [4,18]. 24 x 3600 x G 360n . 
Ho = C (1 + 0.033 cos 365 ) (cos(@) cos(ô) sin(ws) 
Plant Location Radiation kWh/m? day d 
Plataforma Solar Almeria Almeria, Spain 4.82 FWs sin() sin(d)) (2) 
SEGS , California, USA 5.86 Here, Gsc is the solar constant, which is the energy of the sun 
Abengoa ISCCS project Ain-Ben-Mathar, Moroco 4.84 bk ied : p P 
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Fig. 4. Average monthly global irradiance in locations selected [27]. 
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distance, outside of the atmosphere. A value for Gsc of 
1367 [W /m?] is used in this paper. Also n is the n-day of the year, 
¢ is the latitude in radians, ô is the declination angle in radians, 
and w, is the sunset hour angle in radians. The declination angle is 
defined by the equation of Cooper [34] as 


5 =23.45 sin (360 (Se) (3) 


365 
And the sunset hour angle, when the incidence angle is 90°, as 
is needed for CSP plants [29], is defined as 
cos(ws) = —tan(@) tan(6) (4) 
Then, it’s assigned a special distribution to the frequency of 
days with a value of the clearness index Kr. The fraction f of the 


days in the month that are less clear than K; is defined with the 
Bendt correlation as function of Kr [29] 


e0KT.min) —e0Kr) 


fK= eOKT,min) —e0KT,max) 6) 
where y is a dimensionless parameter given by 

1.184č—27.182 -e158 
y = —1.498+ 6 
! KT,max—KT,min ( ) 
where é is also a dimensionless parameter given by 
é s KT, max—KT,min (7) 


K T,max —Kr 


The minimum and maximum values of Kr, Krmin, and Krmax, 
respectively, are given by 


Kr, max = 0.6313 + 0.267K;—11.9(Kr—0.75)® (8) 


Krmin = 0.05 ios 

To obtain a daily clearness index, Knight [9] defines daily Ky as 

In Ie = ae) e0 Krmin) 4. (ae) erm] 
Y 


where nd; is the day of the month and ndm is the number of days 
of the month. 

With all mentioned equations solved, artificial months with 
artificial daily radiations (H) are created, where months are 
arranged from lowest to highest H. Although this procedure 
neglects the autoregressive nature of solar radiation, according 
to Larrain et al. [4], it can be used as a first approximation that, 
nevertheless, constitutes an improvement from monthly mean 
computations. 


Kr= (10) 


3.3. Estimation of hourly diffuse and beam radiation 


As CSP plants accept only DNI, diffuse radiation is subtracted 
from global radiation to obtain the beam radiation, which is the 
one we are interested. The daily diffuse radiation (Ha) is defined 
by Erbs correlations in [29], where depending on the sunset hour 
angle (w,), the daily total diffuse fraction is defined as 

For w, < 81.4° 


1—0.2727K7 +2.4495K7?—11.951Ky° 
pad = +9.3879K;*, Kr <0.715 (11) 
0.143, Kr>0.715 


And for w, > 81.4° 


F 1+0.2832Kr—2.5557K7? 
F =? +0.8448K7°, Kr< 0.715 (12) 
0.175, Kr>0.715 


With H and Hg calculated for each day, hourly radiation (J) is 
obtained by the ratio of hourly to daily total radiation (r;),which is 


defined as the following equation from Collares—Pereira and Rabl 
[29] as function of the hour angle (w in radians) and the sunset 
hour angle (w,): 


a I _ it cos(w)—cos(ws) 
"=J z4 (a+b cos) (Sey ae a (13) 
With a and b constants given by 
‘ 607 
a=0.409+0.516 sin( w.- $35) (14) 
b=0.6609—0.4767 sin (w- 780) (15) 


The hour angle (w) is adjusted by the Spencer Time equa- 
tion [35] and the difference from each location longitude to the 
longitude of reference for Chile. 

Then, based on [30] and Liu and Jordan [28] assumption that 
IaļHa is the same as Ip/Ho, where lo is the hourly extraterrestrial 
radiation, the hourly diffuse radiation Ig is obtained with the ratio 
of hourly diffuse to daily diffuse radiation rg, which is defined as 


" Iq T cos(w)—cos(ws) 
a Ha 24 \sin(w;)—(nw;/180) cosw, 


Finally, hourly beam radiation I, is calculated by subtracting Ig 
from I. 


(16) 


4. Technology selection 


The solar power plant considered in the model consists in 
three main circuits: the solar field, through which the heat 
transfer fluid (HTF) circulates; the power block (PB), where water 
and steam are circulated, equipped or not with an auxiliary 
boiler; and the thermal energy storage system (TES). 

The solar field consists of a number of parabolic trough 
collector loops connected in parallel to each other [35]. Each loop 
consists of several solar collectors, which are made of highly 
reflective parabolic mirrors with receiver tubes mounted at the 
focal line of the parabolic surface. The solar collector total area 
varies according to each power plant capacity considered for the 
simulation [36]. The total CSP plant area is calculated according to 
a 1:4 proportion between the solar collector area and the total 
land area [21]. 

The HTF is distributed to the loops in the solar field by means 
of insulated pipes. A single-axis tracking mechanism allows the 
collectors to follow the sun from sunrise to sunset [37]. The 
sunlight reflected within the trough is focused along a line 
running the length of the trough. In order to collect this heat, a 
pipe is positioned along the length of the trough at its focus anda 
heat collection fluid is pumped through it [18]. The tube receiver 
is designed to be able to absorb most of the energy focused onto it 
and must be able to withstand the resultant high temperature. 
The solar array of a parabolic trough power plant consists of 
several parallel rows of parabolic reflectors. The heat collecting 
fluid, which is pumped through the pipes along the length of each 
solar trough, is typically synthetic oil, similar to engine oil, 
capable of operating at high temperature. During operation it is 
likely to reach between 300 °C and 400 °C [18]. 

After circulating through the receivers of the solar field, the oil 
is delivered to the power block, consisting of a heat exchanger 
train (preheaters, steam generators, superheaters, and reheaters), 
a steam turbine and electricity generator, a condenser, a wet 
cooling system (supplied from sea water pumped to the CSP 
plants site according to the volume that each power plant 
capacity requires [17]), and other auxiliary equipment [37]. The 
heat exchanger extracts the heat from the synthetic oil to raise 
the steam temperature in a separate sealed system. The steam is 
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then used in a Rankin cycle to drive a steam turbine generator to 
produce electricity. The heat collecting fluid is then cycled back 
through the solar collector field to get in more heat. 

Considering solar power alternatives, the parabolic trough 
collector (PTC) power plant is seen as the most convenient for 
solar electricity generation due to its better performance than 
other alternatives and its status as the most developed technology 
[11,12]. The high quality and quantity of specific available data 
about PTC components and processes are the reasons why this 
technology was selected in the present study. One of the previous 
researches [4], in which this study is based in, selected the 
parabolic trough technology with a direct steam generation 
(DSG) process in order to decrease costs related with heat 
exchangers. Even though, we consider the traditional process 
using synthetic oil since the heat collection fluid, as described 
above, is a more mature and developed technology [18]. 

Since the solar radiation displays daily, seasonally and 
multi-annual variability, the configuration of the CSP plant 
counts with a storage system that consists of two tanks of 
molten salts sized for working 7h at full load without solar 
irradiation [38]. At the design-point conditions, collectors 
receive the solar irradiation while they are cooled by means 
of the HTF. Thus, solar collectors produce the nominal mass 
flow of HTF at its nominal temperature. Because the solar field 
is oversized according to the solar multiple, part of the HTF 
mass flow is sent to the steam generator. The surplus mass flow 
is directed to the thermal storage system. In the charging mode, 
the molten salts are pumped from the cold tank, pass through a 
heat exchanger where they are heated and finally go to the hot 
tank at a slightly lower temperature than the nominal HTF 
temperature. When the solar irradiation is lower, i.e. at the 
sunset, production of HTF decreases. If the steam turbine is 
required to work at full load, the thermal storage system begins 
discharging the energy stored in the hot tank to the HTF. 
Therefore, part of the HTF by-passes the solar collectors and 
is warmed up by the molten salts. Molten salts are pumped 
from the hot tank; pass through the heat exchanger where they 
warm up the HTF and finally are discharged into the cold tank. 
Both flows, coming from the solar collectors and from the 
molten salt heat exchanger are mixed and then directed to 
the steam generator [18]. 

Also, in one of the plant arrangements, an additional gas- 
fired boiler will be incorporated to the system in order to 
supply backup heat during night, or stabilize steam properties 
to maintain a constant power output for cloudy days whenever 
the TES is not sufficient [39]. As the fossil fuel backup required 
for maintaining constant generation is function of the local 
conditions (radiation levels, clarity, temperature, and TES load 
level), the plant performance can significantly vary for different 
locations [4]. 


5. Levelized energy costs modeling 


The core of this study is to compare the LECs given by different 
sizes of solar power plants placed in distinct locations of northern 
Chile. The LEC represents the generated electricity costs that 
include the initial capital, return of investment, and variable 
costs. It consists in the division of the total costs by the total 
energy produced. The result is referred as the levelized cost of 
energy, which is widely used to compare competing energy 
sources and is normally calculated using U.S. dollars. The formula 
for calculating LEC is the following: 


AIC 


Eannual 


LEC = (17) 


Table 2 
Values of economics parameters used in the simulation. 


Parameters Value 


Solar field and block of power investment [40] 6,059,461.70 US$/MW 
Storage cost [17] 1,110,425.05 US$/MW 


Backup investment cost [21] 134.960 US$/MW 
Annual O&M cost [40] 70 US$/kW yr 
Price for gas [14] 0.0164 US$/kWh 


Discount rate (r) 10% 

Life time of the plant (n) 25 yrs 

Land rental costs [41] 305.3 US$/ha month 
Annual costs for piping and pumping Vertical: US$ 0.068/ 


100 MM m?/yr [17 m? x 100m 
y 
Horizontal: US$ 0.080/ 
m? x 100 Km 
Table 3 
Values of technical specifications used in the 
simulation. 
Description Value 
Solar field efficiency [42] 0.492 
Thermal energy storage efficiency [42] 0.996 
Power block efficiency [42] 0.4 
Pumping consumption [21] 10% 


where Egnnua is the annual electricity generation and AIC is the 
annualized cost, which is calculated as follows: 


(1+r)"r 
Ca +nr)"-1 


where r is the discount rate, n is the life time of the plant in years, 
and TICVC corresponds to the total investment costs and variable 
costs. In the total investments are considered the costs of the 
solar field, block of power, energy storage system, the gas-fired 
backup system, pipes and bombs for extracting sea water, and for 
cooling processes (annual costs for piping and pumping vary 
according to the distance from the plant location to the sea). 
Variable costs include annual operation and maintenance cost, 
price of natural gas used in the backup system, land rental cost 
and annual costs for pumping energy consumption, and the last 
one is specified in Table 3. 
The variables used for the analysis are resumed in Table 2. 


AIC =TICV (18) 


6. Simulation 


The simulation algorithm is written and developed in MatLab 
commercial software. The inputs used for the analysis are: the 
estimation of daily radiation for each selected site and the 
investment and variable costs previously explained. As output, 
the simulation provides: the values of LEC, the proportion of solar 
energy from the total energy generated, and the capacity factor 
for every plant size and location. The power plant capacities 
analyzed vary between 5 and 200 MW, increasing 5 by 5 MW 
from the minimum value until the maximum value stipulated, 
which give us a number of 40 distinct power plants capacities 
being evaluated by the simulation model, for each location. To 
decide about the correct power capacity range to analyze, it was 
necessary to gather information of other plants currently in 
operation [21]. The optimal LEC in those plants was found to be 
around an installed capacity of 150 MW. Therefore, it was used a 
range between 5 MW and 200 MW for the present analysis. 

The algorithm uses the geographical and meteorological data and 
the solar field area specifications to calculate the corresponding 
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solar time and angle of incidence of the solar radiation over a 
collector operating on sun-tracking mode. These results, together 
with calculations of the solar field efficiency, are used to obtain the 
useful thermal power energy delivered by the solar field. Depending 
of the thermal energy delivered, the decision of the plant’s operating 
mode is determined (TES system and/or backup system operation). 
Consecutively, the corresponding thermal power is sent to the 
power block, which is calculated according to the correspondent 
specifications (energy storage efficiency and/or auxiliary boilers 
efficiency). These results pass through the efficiency specifications 
of the power block (PB) and then, are used to calculate the gross 
electric power generated, the plant’s parasitic consumption and, 
finally, the net electricity production. 

The code is programmed in such a way that the main function 
calculates, for every location and power plant capacity, the energy 
production of every hour of a single day. After a given hour is 
calculated, certain results (like the energy stored in the TES 
system) are saved as input for the calculation of the following 
hour. The finality is to calculate the final energy production of an 
entire year. Also, the algorithm separates the year into a summer 
and a winter seasons. This classification allows to differentiate the 
energy demand peak hours of these two periods. In summer, 
the peak hours are reached from 19:00 to 23:00. During winter, 
the peak hours occur from 18:00 to 23:00. The algorithm does not 
distinguish weekends and holidays, therefore the peak hours are 
applied for every day of the year according to its correspondent 
season. 

An important aspect of the simulation model developed is that 
the minimum level of production of the PB can be set by the user 
[21]. In this work, the minimum load fraction (ML%) is defined at 
30% of the nominal capacity of each plant size. This minimal 
operation level will determine the decision of when the energy 
storage or the auxiliary boiler will work and how much energy 
will be produced by them. When the radiation level is insufficient 
to produce at least at the minimal level stipulated, the thermal 
storage and/or the auxiliary power system are employed in 
parallel to produce the supplement of energy required to reach 
the ML%. However, in peak hours, the plant is considered to 
produce at nominal capacity, which means that is not necessary 
to the solar field to produce lower than the ML% to turn on the 
storage or backup supply. 

The following steps describe the most important aspects of the 
algorithm used for calculating the annual energy production 
according to the ML%: 


(i) The direct solar radiation is calculated in order to determine 
the energy flux received by the solar field. 

(ii) If the thermal flux is higher than the required value to 
activate the plant at nominal capacity, the plant produces 
electricity just with solar energy. The excess energy is stored 
and the auxiliary system is off. 

(iii) If it is not a peak hour and the energy production is not 

enough to reach 100% of the plant capacity, but it is higher 

than the minimum level to work at ML% of the capacity, the 
plant works with solar input from collectors. The stored 
energy does not change, and the auxiliary system is off. 

If the collector steam production does not reach the mini- 

mum level of ML% in hour of regular consumption, there are 

two possibilities: 

(a) The energy storage is able to complement the collector 
field to reach the steam mass production necessary to 
activate the plant at the ML%. The plant works at ML% of 
its capacity and the energy in the storage tanks is 
consequently reduced. 

(b) The stored energy in the tanks is not sufficient to reach 
ML% of the nominal capacity of the plant. In this case, the 


(iv 


Dan 


auxiliary system is enabled to complement the produc- 
tion of the collector field and the tanks to activate the 
plant at a level of MLZ. 
(v) If it is a peak hour and the energy production is not enough to 
reach 100% of the plant capacity, a similar procedure as in (iv) 
will be performed in order to reach the full load capacity. 


In addition to the previously described points, the algorithm 
searches and reserves the 10 consecutive days with lowest 
radiation levels in the year for performing maintenance. This 
procedure allows minimizing the costs of turning off the opera- 
tion of the entire solar plant. 

Also, a second simulation modality is performed, where no gas 
auxiliary system is used. For this to be done, the algorithm to 
produce energy is changed. The previous steps remain the same 
until step (iii), and step (iv) is changed: if the collector steam 
production and the energy storage do not reach the minimum 
level of MLZ, then no energy is produced, the plant is off. Step (v) also 
changes in this second simulation, no peak hours are considered. 


7. Results 


For the 16 locations, the LECs evaluated for the first plant 
configuration (equipped with TES and back-up system) vary 
between 18.72 ¢US$/kWh and 21.85 ¢US$/kWh. The three best 
locations, according to the lowest LECs criteria, are Pica, Coya Sur, 
and Calama, as shown in Table 4. These LEC results are consistent 
with several references in the literature for solar thermal energy 
prices, such as [14,13,15]. 

Fig. 5 shows the results for hourly beam radiation (J,) in 
October 22 for Pica. As shown, radiation peaks to levels of 
3.68 MJ/m? h. However, this radiation levels are not the highest 
among all locations evaluated. This situation suggests that the 
consistency in time of the high levels of radiation allow these 
three locations to operate within a less expansive range, in 
opposite of other ones with higher radiation peaks, but no 
consistency in time of solar incidence levels. 

Energy results in Table 4 show that Calama has the highest solar 
fraction of the annual electricity generation between the three 
locations, while Coya Sur uses the lowest fraction. Despite Calama 


Table 4 
Simulation results: LEC, solar production percentage and capacity factor percentage 
values of top three locations, Pica, Coya Sur, and Calama. 


Parameter Pica Coya Sur Calama 
LEC [ ¢US$/kWh] 18.72 18.82 18.91 
Solar production [%] 77.5 74.83 82.54 
Capacity factor [%] 58.79 58.07 60.04 
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Fig. 5. Hourly radiation beam for Pica. 
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has a higher solar fraction than Pica and Coya Sur and a higher 
capacity factor than both others, its LEC is higher than the others. 
This is explained mainly by Calama’s distance with the sea and 
altitude, what increases piping and pumping investment costs. 

Results also show that a higher nominal power capacity is 
normally related to a lower LEC, as can be observed from a 
200 MW plant in Pica, which have the lowest levelized costs. The 
same power capacity configuration was found for Coya Sur and 
Calama’s best LEC values. 

The results without auxiliary system show that the three best 
locations for solar energy in Chile are: Calama, with a LEC of 
27.55 ¢US$/kWh and capacity factor of 39.57%; San Pedro de 
Atacama, with a LEC of 28.04 ¢US$/kWh and a capacity factor of 
39.25%; and El Tatio, with a LEC of 28.29 ¢US$/kWh and a 
capacity factor of 41.85%. All three achieve these results with 
200 MW of installed power capacity. 

According to [4], the methodology used here for obtaining 
hourly beam radiations has a significant difference compared 
with another methodology that uses identical days with monthly 
mean radiation. This occurs because the second alternative has a 
significant loss in accuracy, which can impact on the usability of 
the collectors. Therefore, the results previously mentioned count 
with more accurate values of available radiation. 


8. Conclusions 


Considering global and local concerns and objectives regarding 
the environment and fossil fuels depletion, finding substitutes to 
conventional energy sources, such as solar energy has become a 
global issue. 

It is known that the SING region, in northern Chile, will be 
highly affected by the increase in energy demand due to mining 
and industrial activity, that are very developed in this area. Solar 
resource in the Atacama Desert becomes, therefore, in an attrac- 
tive alternative for energy generation in the north of the country. 
Considering this scenario, the present study presents an analysis 
of levelized energy cost (LEC) for different power capacities of CSP 
plants placed in distinct locations in northern Chile. 

The results show that CSP plants can be implemented in 
Atacama Desert. For obtaining the optimal LEC, the recommended 
installed power capacity is 200 MW, where the location of Pica 
has the lowest LEC (18.72 ¢US$/kWh). Despite its large distance 
to the seacoast and radiations levels there are significantly high, 
allowing to obtain attractive levelized energy costs. LEC results 
are consistent with several references in the literature for solar 
thermal energy prices, such as [14,13,15]. 

The simulation of a plant functioning without auxiliary system 
shows that levelized energy costs are significantly higher than for 
plants that include this system, with differences of at least 
9 ¢US$/kWh. 

Due to shortage of both ground and surface water in Atacama 
Desert, sea water consumption was considered for feeding the 
cooling system in the simulations. It is noticed that this feature 
lead to significant cost surges in the analysis. 

It is important to remark that the methodology used for 
estimating the available radiation has several limitations. First, 
the autoregressive nature of solar radiation was considered only 
indirectly, since the method is based on empirical correlations 
formulated by observing large data set. Therefore, they might 
include autoregressive effects. Second, the present work intends 
to model the conditions of northern Chile, where high clearness 
index results in a large number of clear days in any given month. 
As a result, the model is not valid for sites that display from 
moderate to large number of cloudy days in a given year. Finally, 
our procedure uses monthly averages in order to develop an 


artificial set of days for each month, and then uses daily totals to 
create an hourly database by means of a clear day model. The 
results could improve in accuracy if a database of hourly radiation 
empirically measured were currently available. 

The authors are currently working in researches that analyze 
different methodologies to reduce costs of CSP plants in Chile, 
such as the trade of carbon emission reduction units. 
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